Controllable Coupling between Flux Qubits 
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We propose an experimentally realizable method to control the coupling between two flux qubits. In our 
proposal, the bias fluxes are always fixed for these two inductively-coupled qubits. The detuning of these two 
qubits can be initially chosen to be sufficiently large, so that their initial interbit coupling is almost negligible. 
When a time-dependent magnetic flux (TDMF) is applied to one of the qubits, a well-chosen frequency of the 
TDMF can be used to compensate the initial detuning and to couple two qubits. This proposed method avoids 
fast changes of either qubit frequencies or the amplitudes of the bias magnetic fluxes through the qubit loops, 
and also offers a remarkable way to implement any logic gate as well as tomographically measure flux qubit 
states. 

PACS numbers: 03.67.Lx, 85.25.Cp, 74.50.+r 



Introduction. — Superconducting Josephson junction cir- 
cuits currently provide one of the best qubit candidates, and 
experiments have been performed for charge, flux, phase, 
and charge-flux |1) qubits. Quantum coherent oscillations 
and conditional gate operations have been demonstrated us- 
ing two-coupled superconducting charge qubits ||2j]. Further, 
entangled macroscopic quantum states have been experimen- 
tally verified in systems of coupled flux [3], and phase 0,|5| 
qubits. 

Quantum computing requires that the interaction between 
different qubits can be selectively switched on and off. This 
is an extremely difficult and important issue. Several schemes 
have been proposed to realize controllable couplings and local 
qubit operations. One is a controllable couplingby dynami- 
cally tuning the qubit frequencies, e.g., in Refs. |£||£|lLI§LE9I- 
This tunable proposal requires that different qubits have the 
same frequencies (i.e., resonant interaction) when they are 
coupled. When they are decoupled, one of their frequen- 
cies should be suddenly changed by an external bias vari- 
able such that two coupled subsystems have a larger detun- 
ing (i.e., non-resonant interaction). The second approach uses 
switchable couplings in charge qubit circuits by changing the 
bias magnetic flux, e.g., in Refs. ifToL HHl . In practice, the 
switching time of the magnetic flux should be less than the in- 
verse single-qubit Josephson energy (less than a nanosecond), 
which is a challenge for present experiments. The third pro- 
posal requires additional subcircuits, e.g., in Refs. fL?! fiUl . 
These additional elements increase the complexity of the cir- 
cuits and add additional uncontrollable noise. 

To easily switch on and off the coupling among qubits is 
one of the most important open problems in quantum informa- 
tion hardware. Here, we propose a way to overcome this se- 
vere problem plaguing experiments. Specifically, we present 
a proposal on how to achieve a controllable interaction be- 
tween flux qubits by virtue of time-dependent magnetic fluxes 
(TDMFs). Here, we make the same assumption as in the de- 
coupling experiments JH0.0], which require the two qubits to 



be in the large detuning regime. However, here, the two-qubit 
coupling and decoupling are controlled by the frequency (not 
the dc component) of the applied TDMF. So we completely 
avoid having to quickly change the bias magnetic flux — a se- 
vere problem faced by many previous proposals for supercon- 
ducting qubits. Moreover, our proposal does not require addi- 
tional subcircuits. These merits also make our proposal poten- 
tially useful for a variety of other types of qubit experiments, 
and could solve a central problem in this field. 

Controllable Hamiltonian. — Two flux-quibts interact with 
each other through a mutual inductance M, as shown in Fig.[0 
Each qubit loop contains three junctions, one of them has an 
area a times smaller than that of the two identical junctions. 
The larger junction in the Zth qubit loop has Josephson en- 
ergy Ejj (I = 1, 2). The gauge-invariant phases (of the 
two identical junctions and the smaller one) in the Ith qubit 
are tp^', <pip, and <pip. We assume that a static (dc) mag- 
netic flux $e and a time-dependent magnetic flux (TDMF) 
$>e\t) = Ai cos(ujcH) are applied through the Zth qubit. Us- 
ing the phase constraint condition through the Zth qubit loop 

+ (27r$?V$o) + (27r#(*)/$o) = 0, the total 
Hamiltonian of the two qubits can be written as 1 14] 

2 2 

H = J2(Hi + H^)+ H lrn +H c + H A . (1) 

1 = 1 l=^m=\ 

Here, the single qubit Hamiltonian is Hi = Pp,/2Mp t i + 
P^/2M Qtl + 2E u (l - coB^cospj?) + aBj,i[l - 
cos(2y>p? + 27r/|)] with the redefined phases (fp = ((p± + 
y>2^)/2, <Pq = (ipi — <p-p)/2, and reduced bias mag- 
netic flux fi = <&e /$o- The effective masses are Mq.i = 
2($ /27r) 2 Cj ! i and M P>1 = (1 + 2a)M QtU which corre- 
spond to the effective momenta Pqj — —ifod/dtfQ and 
Pp i = —ihd/d(p^p. The capacitances in the Zth qubit loop 
satisfy the condition c[ l) = C { 2 1) = C u and C { 3 1) = aC. u . 
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The Hamiltonian i^ = -(A t /2) (i® +i/3P P ^ e~ iu « h + 
H.c. represents the interaction between the Zth qubit and its 
TDMF. Here, the parameter = 27rat4' ) /[*o(l + 2a)], and 
J3 = — (27rai?j j ;/$o) sin(2( / jp' 1 + 2-nfi) is the supercurrent 
through the smaller junction of the Zth qubit without applying 
the TDMF. So a TDMF-controlled single-qubit rotation can be 
realized by the Hamiltonian Hf) . The qubit-qubit interaction 
Hi m , controlled by one of the TDMFs (<S>i 1} (t) or $i 2) (*)), 
can be described by H lm = -/^I^e"^ * cos(2<pp + 
2-Kfi) + H.c, where = M(27r/$ ) 2 (^C'i£j,i/2C , j,0, 
and = C m £<(4rV C i W ' ) ) sin ^| m) is the qubit loop- 
current of the mth qubit without an applied TDMF, and 
C^ 1 = Y^=i{^i)~ l - However, the qubit-qubit interac- 
tion He, controlled by simultaneously applying two TDMFs 
(<f>^\t) and $i 2 \t)) through the two qubits, respectively, 
is H c = Bltf =1 <f> { e l) (t)cos(2ip { p l) + 27r/i) with B = 
M {2tt/<S> ) 4 (CiC 2 Ej A E^ 2 /Cj A C^ 2 )- The Hamiltonian 
Ha — MI^'I^ 2 ' denotes an always-on interaction between 
the two flux qubits, without applying the TDMF. 

In the two-qubit basis {|ei), |<?i)} <£> {\g 2 ), |<72)}> where 
\gi) and |e/) are the two lowest eigenstates (ground and first 
excited states) of Hi, Eq. Q can become 11511 

i=i i=i 
E (n{^f^ ) +H.c)(e ta ? )t + e- 1 "? ) *) 

l=£m=l 

- E (^^Vf)e-^* + H.c.) 

+ (A^JV^+A^JVf +H.c.) . (2) 

Here, the terms /t* ct^ W* and flfjjj ct^ct^W'*, as 
well as their complex conjugates, have been neglected by 
considering energy conservation. The operators of the Zth 
qubit are defined as a) = |e/)(e ; | - o+ = 

|e;)(g/|, and a_ = \gi){e{\. The qubit frequency u>i in 

Eq. (13 can be expressed as uj t = ^21^(<pi l) - $ /2) 2 + tf 

with the loop-current jW and the bias flux $e . Here, 
the parameter denotes the tunnel coupling between two 
wells in the Zth qubit |16|. The controllable coupling 
constants are m = Ai(ei\(I { 3 l) + ipP P ,i)\gi)/2, fi£J = 
Atpi(ei,g m \^coB{2(pf + 27r/,)|fli,c m )/2, and fig] = 
A i A(e i ,e m |/( m ) C0S (2^ ) +27r/0|ft,5m)/2. The hard-to- 
control parameters are Ai = M{e\,g 2 \I^ I^\g\,e 2 ), and 
A 2 = M(ei,e 2 |/ (1) / (2) |gi,g 2 )- It is not difficult to derive 
that k; = = nil] = when no TDMFs. Then, since 

b an im 

both bias magnetic fluxes fi are near 1 /2 (the optimal point 
is at /; = 1/2), the Hamiltonian (|3 can revert to the case in 
Refs. 1 3, 6], where the Pauli operators are defined by the states 
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FIG. 1 : (Color online) Two superconducting flux qubits are coupled 
through their mutual inductance M. Each qubit loop includes three 
junctions. The tunable interaction between two qubits can be realized 
by changing the frequency of the external magnetic flux <J>e (t) (Z = 
1, 2) through the Zth qubit. 



of the two potential wells. 

Decoupling mechanism and logic gates. — We assume that 
the two qubits work at the fixed frequencies u>\ and ui 2 , which 
mean that their reduced bias magnetic fluxes // (Z = 1, 2) and 
the frequency difference A = u)\ — oj 2 remain fixed. If the de- 
tuning A is initially chosen to be sufficiently large (such that 
it satisfies the condition: |A| > \\i\/h = \X 2 \/h = \X\/H), 
then the two qubits can be approximately treated as two de- 
coupled subsystems |7] when the TDMFs are not applied. 

By applying the TDMF with the frequency-matching con- 
dition u>c — 0Ji, we can easily derive from Eq. that any 
single-qubit operation of the Zth qubit can be performed via 
the dynamical evolution ui (6i,<j)i) = exp[i#;(e _ ^ , o-f / ' 1 + 
H.c.)]. Here, 6*/ = \ki\t/H depends on the Rabi frequency 
\ki\/H and duration r; <f>i is related to the TDMF phase ap- 
plied to the Zth qubit. For example, 7r/2 rotations of the Zth 
qubit around the xory axes can be performed by Z7 C (9i ,<t>i), 
with t = Htt/A\ki \ and 4>i = tt or 7r/2. Unless specified oth- 
erwise, hereafter, we work in the interaction picture, and all 
non-resonant terms have been neglected because their contri- 
butions to the transitions between different states are negligi- 
bly small 0. 

To couple two qubits with the assistance of the TDMF: i) 
a TDMF needs to be applied through one of the qubits, and 
its frequency should be equal to the detuning (or summation) 
of the two qubit frequencies; ii) the reduced bias flux 111811 on 
the qubit, which is addressed by the TDMF, should be near but 
not equal to 1 /2. Without loss of generality, below, the TDMF 
is assumed to be always applied through the first qubit, so the 
bias of the first qubit is fx = 1/2 ± e, with small e; however, 
the bias for the second qubit is taken as f 2 = 1 /2. 

Considering two new frequency-matching conditions in 
Eq. (|2jl, produces two different kinds of Hamiltonians for 
implementing two-qubit operations with the assistance of 
TDMF. One is H x = 0$ cr W + H.c, with the condition: 
lu 1 —lu 2 ± uic = 0. Here, the sign is positive (negative) when 
A < (A > 0). Another one is H 2 = o { + ] <t { + ] + H.c, 

when the frequencies satisfy the condition: u>i + u> 2 — u>c — 
0. Using the Hamiltonian Hi and H 2 , two qubit gates can 
be implemented. For example, a TDMF is applied through 
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the first qubit with its frequency w' 1 ', satisfying the condition 
lui — LU2 — wi = (hereafter, we assume u)\ > UJ2 without 
loss of generality). The Hamiltonians Hi can be transformed 
to Hi = H XY = -(\n[ 1 2 ) \/2){a x 1) <j x 2) + 4 1} 4 2) l when 
the phase of is set to it by an applied TDMF. Then an 
ISWAP gate [jjj, denoted by C/ iS , can be realized by H X y 
with an evolution time t = irh/ (2\^ 2 ] \). The CNOT gate 
can be constructed by combining the ISWAP gate with a few 
single-qubit operations. 

Experimentally, it is found that the always-on coupling 
strength \\\/h is about several hundred MHz, e.g. \X\/h ~ 
0.4 GHz in Refs. 0,0], the detuning A/2tt can be - 1 to 10 
GHz. This means that the ratio |A|//iA cannot be infinitesi- 
mally small, and the always-on interaction needs to be con- 
sidered when all TDMFs are switched off. Up to the first or- 
der in |A|/fiA, the effect of the always-on interaction, with- 
out the TDMF, can be described by an effective Hamiltonian 
H E = (|A| 2 /ftA)[| ei )( ei | ® \g 2 )(g 2 \ - \ 9l )( gi \ ® |e 2 )(e 2 |]. 
Here, fast oscillating terms have been neglected. 

Entangled states and tomographic measurement. — Entan- 
gled states can be easily generated in this circuit. For example, 
if the first and second qubits are initially in the excited state 
|ei) and the ground state |g 2 ), then using the Hamiltonian 
the system evolves to \$>i{t)) = cos(|fij^ |i//i)|ei, <7 2 ) — 
ie~ %e sindS!^ e 2) with the phase 9 of fi^ . It is ob- 

vious that the Bell states | ip±) = (\gi, e 2 ) ± \ei, g 2 ))/V2 can 

be generated with ti = 7i7r/(4|f2^|) bv setting the TDMF 
such that 9 = n/2 or 3n/2. Similarly, if both qubits are in 
the ground states \gi) and | (72) , then another two Bell states 
l*±) = (151)52) ± |ei, e 2 ))/"\/2 can also be obtained with 
< 2 = 7wr/(4|Q^)|) b Y setting the phase 9' of Q$ as n/2 or 
3ir/2 through the Hamiltonian H 2 . 

State tomography allows us to experimentally determine a 
quantum state 12011 . Qubit state tomography can be imple- 
mented by measuring the supercurrent through the qubit loop, 
which is inductively coupled to, e.g„ a dc SQUID magne- 
tometer or high-quality tank circuit Q. For the 2th qubit in 
the qubit basis {\gi), |e;)}, its loop-current operator can be 
written lEHl as 

/CO = /« = a ( cT« +bl \ ei )( gi \+ b *\ gi ) {ei \ (3 ) 

with ai = (ei\l"'\ei) and bi = (ei\l"'\gi), when the bias 
fi is near (but not equal to) 1/2. However, at the optimal 
point fi = 1/2, the supercurrent operator in Eq. Q can be 
reduced lEHl to /CO = f$ = bi a x \ with a real number bi 
and the Pauli operator ai = |e;) (gi \ + \gi) (e;|. 

If the simultaneous joint measurement of two qubits can 
be performed in flux qubit circuits as in phase circuits 
then single qubit operations are enough to realize the fif- 
teen different measurements [20] on the two-qubit states 
p = (1/4) J2i j r i,j a i ® CT j 2 '' w ith the Pauli operators af' 
(i, j = x, y, z and I = 1, 2) and the identity operator o\ , 
where roo = 1 by normalization. The loop-current operator 
for the first qubit is given in Eq. by setting I = 1 due to the 



(2) (2) 

assumption /1 7^ 1/2, but it is reduced to I x ' = b2d\ for 
the second qubit with / 2 = 1/2. So the fifteen measurements 
on state p are given as /^ and ^cj 2 ^ (denoted as single-qubit 
measurements), as well as 6 2 /j ® o~j^ (called two-qubit or 
joint measurements), with i, j — x, y, z, Iy = Y^li Yi, 
and ijp = z\lx It is clear that three measurements 
(Ix'K b2cr x 2 \ and b2llP ® <J X 2 ^) on the input two-qubit state 
p can be directly performed. Other twelve measurements can 
be equivalently obtained by measuring (I ( x \ b 2 a x 2 \ or both 
of them at the same time) on the rotated state p. For example, 
7r/2 single-qubit rotations Yi around the y axis for the first 
qubit and Z 2 around the z axis for the second qubit are si- 
multaneously performed on the state p, then the measurement 
& 2 /k 1 ' ) <£> o~ x 2 ^ on the rotated state Y\Z%pZ\Yy is equivalent to 
the measurement 6 2 /y 1 ' ) <S> Uy 2 ' on the original state p. Simi- 
larly, other joint measurements can also be obtained. Finally, 
for the fifteen measured results, we solve a set of equations 
for the parameters r\j, and a two-qubit state is reconstructed. 

If only a single-qubit measurement can be made at a time, 
besides six single-qubit measurements mentioned above, a 
suitable nonlocal two-qubit operation 12011 is required to ob- 
tain the coefficients (e.g., r V:Z ) of the nine joint measurements 
on the state p. Here, this is an ISWAP gate Uis, which can 
be implemented as described above. For example, if an oper- 
ation ?7is is made on the input state p, then the loop-current 
of the second qubit should be {I {2) ) = Tr(J7 iS p I {2) ) = 

— 6 2 Tr(p<7y 1 ' ) <g> o^ 2 ') = —b2T y . z , and the coefficient r VjZ is 
determined. Combining the ISWAP gate and single-qubit op- 
erations for two qubits, the coefficients of other eight joint 
measurements can also be determined by only measuring the 
loop-current Ix ■ 

Tomographic ally measured states are different for com- 
pletely decoupled (CD) and large detuning (LD) two-qubit 
systems after two-qubit states are created, if we consider a 
duration t before measuring the generated two-qubit states. 
As an example, a schematic representation of a Bell state 
is given in Fig. |2]for the above two cases. There is 
only the real part for the reconstructed state p = 
in the CD system, shown in Fig. EJa). However, due to 
the effect of He for the LD system, the reconstructed state 
X = e~ lHEt / h \?p + ) (tp + \ e *H E t/h mc i u( jes both real and imagi- 
nary parts, shown in Fig.|2jb) and 0c), respectively. In Fig.|2] 
we consider a longer duration t ~ 10 ~ 9 s; the detuning and 
the coupling constant are, e.g., A ~ 5 GHz, and \X\/h ~ 0.4 
GHz. So if we consider the always-on interaction effect, the 
relative error with these parameters is ~ 0.08 for the non- 
diagonal parts of the reconstructed CD state Here, the 
qubit free-evolution is neglected. In practice, considering un- 
avoidable environmental effects and statistical errors, the ex- 
perimentally measured data should be further optimized by 
other methods |22] . 

Conclusions. — The controllable coupling of two 
inductively-coupled flux qubits can be realized, when 
the large detuning condition is satisfied, by the frequency of 
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FIG. 2: (Color online) A schematic representation of tomograph- 
ically measured state p [with only real part in (a)] for a completely 
decoupled system, as well as a state x [with real part in (b) and imag- 
inary part in (c)] for a two-qubit system with large detuning — after 
considering a duration t before measuring the created state \ip+). 

the TDMF matching/mismatching to the detuning (or sum- 
mation) of the two qubit frequencies; not by changing qubit 
biases (e.g., as in Ref. |7|). Our proposal is also different from 
the coupling/decoupling method by using dressed states 12311 . 
We emphasize that the deviation e from the optimal point 1 /2 
for the reduced bias fi through the first qubit will make the 
decoherence time T 2 short. However, our proposal can work 
for a small deviation, e.g., e ~ 1CP 4 , in which T 2 > 20 ns 
in Refs. 0] or T 2 > 100 ns with spin-echo signals I24ll . At 
this point, the qubit coupling constant fij^ (* — 1> 2) can 
reach 1 18] about several hundred MHz. Based on numerical 
estimates l2lll . the longest time for the single-qubit operation 
Zi is less than 5 ns, so the tomographic measurements can 
be performed within T%. Our proposal can be scalable to a 
chain of many inductively coupled flux qubits, if all of qubits 
satisfy the large detuning condition. We need to note: i) 
we can use one LC circuit as a common information bus to 
couple many qubits, with the qubit-bus coupling controlled 
by externally variable frequencies 1 25 ] ; ii) this circuit can be 
modified to work at the optimal point; iii) this method using 
frequency-controlled couplings can be applied to control 
one-junction flux qubits. It can also be modified to control 
phase, charge-flux, and charge qubits. 
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